NASA Technical Memorandum 105238 
ICOMP-91- 16; CMOTT-91-06 


///-3y 




A k-e Modeling of Near 
Wall Turbulence 


(NA5A-TM lu523u) A k-LPAILON MODELING OF N91-324 

NEAR WALL TURBULENCE (NASA) 14 p CSCL 200 

Unci us 

G 3/34 0045757 


Z. Yang and T.H. Shih 

Institute for Computational Mechanics in Propulsion 
and Center for Modeling of Turbulence and Transition 
Lewis Research Center 
Cleveland, UFio . 


Prepared for the 

Fourth International Symposium on Computational Fluid Dynamics 
Davis, California, September 9-12, 1991 



IWVSA 



A k-e MODELING OF NEAR WALL TURBULENCE 


Z. Yang and T.H. Shih 

Institute for Computational Mechanics in Propulsion and 
Center for Modeling of Turbulence and Transition 
NASA Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 


A Jfe — e model is proposed for turbulent wall bounded flows. In this model, the turbulent 
velocity scale and turbulent time scale are used to define the eddy viscosity. The time scale 
is shown to be bounded from below by the Kolmogorov time scale. The dissipation equation 
is reformulated using this time scale, removing the need to introduce the pseudo-dissipation. 
A damping function is chosen such that the shear stress satisfies the near wall asymptotic 
behavior. The model constants used are the same as the model const ants in the commonly 
used high turbulent Reynolds number k-e model. Thus, when it is far away from the wall, 
the proposed model reduces to the standard k-e model. Fully developed turbulent channel 
flows and turbulent boundary layer flows over a flat plate at various Reynolds numbers are 
used to validate the model. The model predictions are found to be in good agreement with 
the DNS data. 

INTRODUCTION 

Because of a wide range of scales involved in a turbulent flow, DNS (direct numer- 
ical simulation) is limited to flows of moderate Reynolds number and simple geometry. 
Turbulence modeling is the only viable approach for the calculation of turbulent flows of 
engineering interest. In turbulence modeling, k-e model is the most used model in engineer- 
ing calculations and model constants have been computationally optimized, giving what is 
called the Standard Model (e.g., Spalding and Launder [1], Rodi [2]). The Standard Model 
is devised for high turbulent Reynolds number flows and is traditionally used in conjunction 
with a wall function when it is applied to wall bounded turbulent flows. Since universal 
wall functions do not exist in complicated flow situations, it is necessary to develop a form 
of the k — e model such that the equations can be integrated directly down to the wall. 

Jones and Launder [3] made the first proposal for the low Reynolds number k — e model 
and the equations were successfully integrated to the wall. A number of models have been 
proposed since then. A critical evaluation of the pre-1985 models was made in the paper by 
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Patel et al. [4]. They pointed out that it is important for the damping function in the eddy 
viscosity to satisfy the near wall asymptotics. More recent models could be found in Shih 
[5]. Of the models proposed, a pseudo- dissipation is introduced near the wall. However, 
the definition of the pseudo- dissipation is quite arbitrary, except for the constrains on the 
wall and far away from the wall. In addition, the model constants used are different from 
the Standard Model, making the near wall models less capable of handling flows containing 
both high Reynolds number turbulence and near wall turbulence. 

In the present paper, we try to provide remedies for these two shortcomings stated 
above. In section 2, we list the general form of the low Reynolds number k - e models. We 
then propose in section 3 a new model. The proposed model is validated in section 4 for 
fully developed turbulent channel flows and turbulent boundary layer flows over a flat plate 
at different Reynolds numbers. Section 5 presents conclusions and possible future works. 

EXISTING NEAR WALL k - e MODELS 

In turbulence modeling, the instantaneous quantities of an incompressible flow are de- 
composed into the mean and the fluctuations, i.e. = Ui + u», p = P + p. The mean 
field U{ satisfies the following continuity equation and the Reynolds averaged Navier-Stokes 
equation: 

Ui,i = 0 ( 1 ) 

U{ + UjUij = — — Pi + vUijj— < UiUj > j ( 2 ) 

where the Reynolds stress term, — < U{Uj >, has to be modeled. 

In an eddy viscosity model, the Reynolds stress is related to the mean field by 

2 

- < UiUj >= v T (U itj + Uj t i) - kSij , (3) 

where up is the eddy viscosity and k is the turbulent kinetic energy. From a dimensional 
reasoning, the eddy viscosity is given by 

vt ~ Uthi 

where u t and l t are the turbulent velocity scale and turbulent length scale respectively. 

In the framework of the k — € model, l t ~ k 3 / 2 /e and u t ~ fc 1 / 2 . In the case of near 
wall turbulence, the above is modified by introducing a damping function and a pseudo- 
dissipations €. The general forms for the eddy viscosity and the transport equations for k 
and e are (Patel et al. [4]), 

V T = C pf n~~ (4) 
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( 5 ) 


k + Ujkj = [(u + —)k,j],j~ < UiUj > Uij - e + D 

<Jk 

k + UjCj = [(u + y )€j]j + Cift^PTUijUij - C 2 f 2 j + E (6) 

where 

€ = €-D 

Various models use different sets of model constants c^, <7fc, <r 0 c*, c 2 and different 
damping functions / 1 , / 2 as well as different forms of Z) and Z7. Far away from the wall, 
the damping functions / 1? / 2 approach to 1 and the near wall terms Z?, E are negligibly 
small so that the model equations become the high Reynolds number forms of k — c model. 

THE PROPOSED MODEL 


The turbulent length scale is characterized by the size of the energy containing eddies. 
Near the wall, these eddies would have a size of 0(y). Following Patel et al [4], the turbulent 
velocity field has the following expansions near the wall: 

u 1 - u x y + u 2 y 2 + ... 
v 1 = v 2 y 2 + ... 

w 1 — w x y + w 2 y 2 + ... (7) 

where non-zero in general. Thus, as the wadi is approached, both the turbu- 

lent length scale and the turbulent velocity scale approach zero. However, the turbulent 
time scale, which is given by the ratio of the scale of the energy containing eddies to the 
turbulent velocity scale, approaches to a non-zero value. We expect this time scale to be the 
Kolmogorov time scale because viscous dissipation is dominating near the wall. We thus 
have, for the turbulent time scale, 

T t = * + T fc (8) 

where 
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v 



and the eddy viscosity is 


V T = c^f^kTt ( 11 ) 

0, since we are 

( 12 ) 

(13) 

(14) 

Since T t is always a positive nonzero quantity, the above formulation removes the singu- 
larity of the high Reynolds number k — e model when the wall is approached. / 2 and E are 
of the form used by Jones and Launder [3] for the low Reynolds number model, in which f 2 
is used to capture the final decay of the homogeneous turbulence and E is used to model 
turbulence in the buffer layer. 

Near the wall, the shear stress - < uv > should behave as 0(y 3 ). Thus we would 
require a damping function which has a near wall behavior of O(y). The damping function 
is assumed to be of the following form 

fn = 1-0 - exp(a l2 / + + a 2 (y + ) 2 + a 3 (y + ) 3 + a 4 (y + ) 4 ) (15) 

with the constants are determined by calibrating with the DNS data for 2D channel flow at 
Re r = 180. The calibration gives 

ai = 4xl0~ 3 , a 2 = 5xl0“ 5 , a 3 = -2xl0“ 6 , a 4 = 8xl0“ 8 

The model constants used are the same as those in the Standard Model, i.e. 

= 0.09, Ci = 1.44, C 2 = 1.92, = 1.0, a t = 1.3 

Thus, far away from the wall, / 2 , become 1, the Kolmogorov time scale is much smaller 
than fc/e, and B is much smaller than the other terms on the right hand of e equation, our 
model would reduce to the standard k — e model, ensuring the good performance of our 
model in the high turbulent Reynolds number limit. 

MODEL VALIDATION 

We use 2D fully developed channel flows and turbulent boundary layer flows to validate 
the proposed model. These flows are attractive for model testing, because they have self- 
similar solutions so that the initial conditions do not have to be accurately specified. These 


The transport equation for k remains the same as eq (5) (with D = 
solving the real dissipation.) The transport equation for e is 

e + Ujej = [{v + + {CxurUi^Uij - C 2 f 2 e)± + E 

a € It 


h = 1-0- 0.22 exp(— 
E = vvTUijkUijk 
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flows are very simple, and solutions can be found very efficiently; yet, the effects of the wall 
on turbulent shear flow are still present. In addition, DNS data providing detailed flow 
information is available for these flows. We will compare the model predictions with the 
DNS data. 

The boundary condition for e on the wall is determined by applying k equation to the 
wall, which gives 


€ W — Vkyy, ( 16 ) 

In this study, the following boundary condition for e, which is mathematically equivalent 
to the above but computationally much more robust, is used. 

dk 1 ! 2 

«• - <»> 

An implicit finite difference scheme is used to solve the momentum equation and the 
transport equations for k and e. The coefficients for the convective terms are lagged one 
step in the marching direction, and the source terms in the k find e equations are linearized 
in such a way that numerical stability is ensured. 

A variable grid spacing is used to resolved the sharp gradient near the wall. The grid 
distribution is controlled by Syi/Syi^i = a. Both a and the total number of the grid, N, are 
varied to ensure the grid independence of the numerical results. The marching step size, 
Sx y is also varied to ensure accuracy. It is found that N = 150, a = 1.03, and Sx — 0.05 are 
sufficient to give a solution with a less than 0.5% error. 

Computations are carried out for 2D fully developed turbulent channel flows at Re r = 
180 and Re T = 395. Fig. 1 shows the mean velocity profile for Re r = 180, together with the 
DNS data of Kim et al [6]. In this figure, and in the following figures, both the dependent 
and the independent variables are normalized by u T and u. The agreement between model 
prediction and DNS data is excellent. This may not be unexpected since the constants in 
the damping function are calibrated for this case. We then show the results for 2D channel 
flow at a different Reynolds number, Re T = 395. Fig. 2 to Fig. 5 shows the profiles of the 
mean velocity, shear stress, turbulent kinetic energy and the dissipation, respectively. The 
predictions of Jones and Launder [3] and Chien [7] are also shown. These predictions are 
compared with the DNS statistics for this case (Mansour [8]). Overall, the proposed model 
gives a better prediction. 

Predictions of the mean velocity, shear stress, turbulent kinetic energy and dissipation 
rate for turbulent boundary layer over a flat plate at Re$ = 1410 are shown in Fig. 6 to 
Fig. 9. The model predictions are compared with the DNS data of Spalart [9], Similar 
computations are made for Re$ — 670, where the DNS data is also available (Spalart 
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[9]). The results for Reg = 670 are not shown here, due to space limitation. The model 
predictions for the skin friction are shown in Fig. 10, together with the DNS data of Spalart 
[9] at Reg = 670 and Reg = 1410. 

In the above computations, c& = 1.0 is used. The constant c* was varied in the range of 
0.5 - 3.0 and the solutions were found to be quite insensitive to cjt in this range. As more 
flow situations are tested, the value of Ck could be optimized by fine tuning. 

CONCLUSION 

We have presented a k — e model for wall bounded flows. We have shown that near the 
wall, the time scale is characterized by the Kolmogorov time scale rather than zero. By using 
this time scale and a turbulent velocity scale, the singularity in the standard dissipation 
equation is removed as the wall is approached and the equation can be integrated down 
to the wall. The use of this time scale also renders pseudo-dissipation unnecessary, which 
contains an arbitrary element in its definition. 

The proposed model uses the same set of model constants as that used in the Standard 
Model, and when it is far away from the wall, the proposed model reduces to the Standard 
Model. Thus, the proposed model would be applicable in both the near wall turbulence 
and the high Reynolds number turbulence. 

The damping function is chosen in such a way that the near wall asymptotics for the 
shear stress is satisfied, a feature which is very important in predicting the near wall prop- 
erties (Patel et al [4]). The form of the damping function is assumed to be a function of 
y + with the constants calibrated with the DNS data for a particular case. This y + depen- 
dence of the damping function must be viewed as a drawback of the model, for y + is not 
clearly defined, or even does not exist for complex flows (comer flows, flows with separa- 
tion, for example.) In addition, y + makes the modeling equation not invariant in coordinate 
transformations. Currently, effort is being made to remove these deficiencies. 
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Figure 7. Turbulent shear strew for flat plate boundary layer, Re$ = 1410. 
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Figure 8. Turbulent kinetic energy for flat plate boundary layer, Re§ = 1410. 
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